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Transport through spin-blockaded quantum dots provides a means for electrical control and de- 
tection of nuclear spin dynamics in the host material. Although such experiments have become 
increasingly popular in recent years, interpretation of their results in terms of the underlying nu- 
clear spin dynamics remains challenging. Here we examine nuclear polarization dynamics within 
a two-polarization model that supports a wide range of nonlinear phenomena. We point out a 
fundamental process in which nuclear spin dynamics can be driven by electron shot noise; fast 
electric current fluctuations generate much slower nuclear polarization dynamics, which in turn af- 
fect electron dynamics via the Overhauser field. The resulting intermittent, extremely slow current 
fluctuations account for a variety of observed phenomena that were not previously understood. 
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The opportunity to study spin coherence and many- 
body dynamics in a controllable solid-state setting has 
inspired a wide range of experiments in a variety of 
materials such as GaAs vertically grown and gate- 
defined structures^, InAs nanowirea^, and ^^C-enrichcd 
carbon nanotubes^. In particular, electron transport 
through spin- blockaded double quantum dots^ consti- 
tutes a purely electrical means of probing and manip- 
ulating the dynamics of nuclear spins. Such experiments 
have revealed complex dynamical phenomena, including 
bistability and hysteresis^i^i^i^, switching^i^i^, slow tran- 
sient build-up of current^, and slow oscillations^!^. 

Despite wide interest in these phenomena and their im- 
portance for quantum information processing, progress in 
understanding them has been slow. While there is little 
doubt that nuclear spins in the host material play a cru- 
cial role, the lack of a direct probe of nuclear spin dynam- 
ics requires their behavior to be inferred from electronic 
transport measurements. To meet this challenge, theo- 
retical modeling must be used to complement analysis of 
relevent features in transport data. 

In previous work on spin dynamics in double quantum 
dots, simple models involving a single dynaical variable 
describing the total nuclear polarization have been used 
to explain the origin of feedback in this systeroiSiiii. Al- 
though such models can succesfully account for feedback- 
driven nonlinear phenomena such as bistability and hys- 
teresis, the range of phenomena which they can describe 
is somewhat limited. Here we expand the phase space of 
the model, and describe nuclear spin dynamics in terms 
of two dynamical variables, and s^, corresponding to 
the independent nuclear polarizations in the left and right 
dots (see also e.g. Ref. Il2h . thereby extending the range 
of phenomena that can be analyzed. Time evolution 
is described by trajectories in a two-dimensional phase 
space (sL,Si?), which can exhibit complex dynamics in- 
cluding non-monotonic behavior, limit-cycles, or spirals, 
as illustrated in FigH] (also see Refs. [H and [13 for addi- 
tional examples of complex phenomena arising from two- 
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FIG. I: Nuclear dynamics and intermittent current fluctu- 
ations driven by shot noise of the spin-blockaded current 
through a double dot. a) Hyperfine spin-exchange with nuclei 
mediates transitions between two-electron triplet and singlet 
states, relieving blockade and producing dynamical nuclear 
polarizations (DNP) sl and sr in the left and right dots, 
b) A phase portrait of DNP trajectories, with a fixed point 
(DNP steady-state) positioned near the main diagonal, where 
current is low due to spin-blockade (see FigO. c) Typical sim- 
ulated current trace showing the effects of steady state DNP 
fiuctuations (pink traces in b and c), see appendix. Switching 
between quiet and noisy regions results from excursions into 
the dark stripe sl ~ sr, marked by the shaded regions in 
panel c). 



polarization dynamics in other contexts). 

Nuclear polarization dynamics in spin-blockaded dots 
is driven by carriers passing through the system^. Each 
electron passing through the dot can produce a spin flip 
of the nuclei due to hyperfine exchange with nuclear spins 
in the host lattice, see FigH^. Dynamic nuclear polariza- 
tion (DNP) arises when the up and down spin-flip rates 
are imbalanced, F4. 7^ F_— li^. Since the spin-flip rates 
in the two dots are in general different, their correspond- 
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ing nuclear polarizations Sl and have different time 
dependence, generating an asymmetry between the dots, 
Sl sr- As we shall see, such asymmetry is dramati- 
cally reflected in the time-dependence of electric current. 

In this work we focus on the effects in nuclear po- 
larization dynamics due to the shot-noise arising from 
the discreteness of carriers passing through the system. 
Electrons are injected into the system one-by-one, with 
random spin orientations. While transiting through the 
dots, each such electron may exchange its spin with the 
nuclear subsystem. Crucially, these stochastic spin-flip 
processes comprise an intrinsic source of broadband noise 
that couples to nuclear dynamics. The intensity of this 
noise, which is proportional to the DC current, remains 
nonzero even when the average rates of up and down 
spin- flips are equal: S oc (r+ -|- r_). The resulting DNP 
fluctuations are relatively slow due to the large number of 
nuclear spins in the dots, N « 10*^, which requires many 
electrons to be transmitted through the system before 
DNP can change substantially. 

Another important aspect of the double dot system 
is the complex relationship between the system's inter- 
nal variables and measurable quantities, i.e. between the 
nuclear polarization and electric current. Due to the res- 
onant energy dependence of transition rates, current is 
sensitive to the alignment of energy levels via a number of 
external and internal variables (gate voltages, magnetic 
field, Overhauser fields in each dot, etc). Changes in the 
hyperfine spin flip rates feed back into DNP dynamics, 
giving rise to a variety of interesting nonlinear phenom- 
ena occurring on long time scales, exemplified in FiglTlD,c. 
Numerical simulations based on this microscopic model, 
which is described in detail below, demonstrate how the 
complex long timescale dynamics arise from the stochas- 
tic nature of electron transport. 

In particular, we find that the high frequency noise 
can drive intermittency in electric current resembling the 
multi-scale switching behavior observed in experiments, 
which will be discussed below. In dynamical systems, 
intermittency refers to the alternation of phases of ap- 
parently regular and chaotic dynamics^^. Such behavior 
arises in many physical systems. For example, fluores- 
cence intermittency, or blinking, is commonly observed in 
the optical repsonse of various nanoscale systems, such as 
large molecules or quantum dots, where it signals com- 
petition between the radiative and non-radiative relax- 
ation pathways^. In our system, a commonly observed 
type of behavior is slow build-up followed by intermittent 
switching between "quiet" and highly fluctuating current 
states, illustrated in FiglT]D,c and Figl5^,b. 

Throughout this paper, simulation results are com- 
pared to data from the measurements described in 
Figure shows typical experimental current traces ob- 
served in the regime of moderate magnetic fleld {B = 
200 mT) and with a gate voltage setting where the elec- 
trostatic energy makes lowest two singlet states, one with 
one electron in each dot and the other with both elec- 
trons in the right dot, nearly degenerate (i.e. near zero 




FIG. 2: Time-dependent current, from experimeniji (a), and 
from simulation (b), showing a long build-up lasting several 
tens of seconds followed by steady-state intermittent fluctu- 
ations (blue traces). Fluctuations can be suppressed by a 
small change of detuning (gate voltage) which moves a fixed 
point away from the sensitive region sl ~ sr (red traces), 
c) Fourier spectra of the experimental and simulated noisy 
traces, evaluated in the marked steady-state regions. The 
spectra display a roughly l/f" dependence with Ooxpcr. ~ 2.5 
and aaimui. ~ 2. 



"detuning"). In this case, these "(1,1)" and "(0,2)" sin- 
glet states are strongly hybridized by the tunnel cou- 
pling between the dots (see Fig|3^). The traces were 
taken after a long waiting period which allowed the sys- 
tem to relax to equilibrium. Current displays dynamics 
on a very long time scale, with a smooth transient "slow 
build-up" period lasting several tens of seconds followed 
by a "steady-state" featuring intermittent large ampli- 
tude fluctuations with a correlation time on the scale of 
seconds (blue trace). The fluctuations can be abruptly 
suppressed by a relatively small change of detuning (red 
trace). Similar behavior was observed during slow sweeps 
of magnetic field (shown in Fig|4]). 

Similar-looking fluctuations were reported by Reilly 
et al. as 'blinking' of the Overhauser fleld measured in 
a double-dot which was repeatedly pulsed through a 
singlet-triplet level-crossing^. There, long timescale noise 
correlations were attributed to nuclear spin diffusion re- 
sulting from the dipole-dipole interaction. In contrast, 
below we describe a mechanism where diffusion of the 
net nuclear polarization is not driven by the conventional 
dipole-dipole mediated spin flips, but rather is driven by 
shot noise in the current passing through the system. 

The rest of the paper is organized as follows. In Sec- 
tion |T] we describe the physical mechanism of shot-noise- 
induced multiscale intermittent fluctuations of current. 
Then in Section [IT] we present the mathematical descrip- 
tion of our model for describing the time dependence of 
nuclear polarization and current in spin-blockaded dou- 
ble quantum dots. In Section Hill we present the results 
of simulations based on the model described in SecHIl 
and compare with experimental data. Finally, our con- 
clusions arc summarized in Section IIVI 
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I. TRANSIENTS AND INTERMITTENCY IN 
THE TWO-POLARIZATION MODEL 

A typical behavior, often seen in the data, is a rela- 
tively slow transient buildup of current after which the 
system enters an intermittent state, characterized by al- 
ternation of quiet and noisy behavior. Here we discuss 
the physics of how such behavior can arise naturally from 
the two-polarization model. The key elements of the 
mechanism arc summarized in schematic form in FigUlD, 
which shows a phase portrait of DNP in the {sltSr) 
plane. In our analysis we assume that, via the hyper- 
fine interaction, dynamics are primarily controlled by two 
variables sl and sr that describe independent nuclear 
polarizations in the two dots. The trajectories shown in 
FiglTJj are obtained by applying the ideas of Refs. HI and 
lid to the case of two coupled polarizations while ignoring 
noise; the fixed points associated with these trajectories 
describe steady-state DNP. Due to blockade of the triplet 
state I To ) , current is low in the gray stripe indicated 
along the main diagonal, sl = sr; away from this line 
the finite polarization gradient As = sl — sr mixes | Tq ) 
with the singlet states [see Figl3] and Eq. ^ below] and 
gives rise to enhanced current, which is then only weakly 
sensitive to DNP. 

Intermittency originates naturally as follows. Due to 
asymmetry of the dots, DNP initially moves from the un- 
polarized state into the region sl ^ sr where current is 
insensitive to polarization (pink curve in FiglT})). This 
corresponds to the quiet build-up period in the current 
trace (see Fig 13]). After approaching the nearly symmet- 
ric fixed point, DNP continues to fiuctuatc locally. Here, 
relatively small fluctuations of polarization take the sys- 
tem in and out of the low-current stripe sl ~ sr, re- 
sulting in large amplitude fluctuations of current and its 
apparent 'switching' between high and low values. 

This behavior arises whenever a DNP fixed point re- 
sides near the sensitive region sl ~ sr, irrespective of 
the details of the dynamics near the fixed point. Using 
a realistic model of the stochastic dynamics of electron 
transport and nuclear polarization, we have generated 
current traces exhibiting both the slow quiet build-up 
and long-time intermittent fluctuations. The relationship 
between the intermittent behavior of current and fluctu- 
ating DNP is illustrated in FigHJ;, where corresponding 
regions of low current and sl ~ sr are marked in gray. 
When parameters such as magnetic flcld or detuning are 
changed such that the DNP fixed point moves away from 
Sl Sr, intermittent fiuctuations of current are abruptly 
suppressed (see red line in Figl^b); this picture is consis- 
tent with experiment (sec Figl2^ and Ref. 0)- 



II. MATHEMATICAL FORMULATION 

Wc now turn to a detailed description of the model 
used to generate FigsU] and [21 The relevant energy lev- 
els arc depicted in Fig 131 The states (1, l)s and (0,2)5, 
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FIG. 3: a) Energy level diagram showing hybridization of 
'S'Tot = levels, see Eq.((T]), and relevant transitions for spin- 
blockaded transport, b) Current as a function of nuclear po- 
larizations Sl and sr on the left and right dots. Current is 
suppressed for sl ~ sr due to the vanishing coupling between 
I To ) and | (1, l)s ). Small roman numerals indicate bands of 
high current due to resonance between blocked and unblocked 
levels, shown in panel a). Primed labels indicate analogous 
configurations with | T± ) reversed, c) Simulated trajectory of 
polarization arcing through the smooth region of the current 
diagram with steady-state near the sensitive region sl ~ sr, 
similar to the schematic in FiglT] Arrows indicate the direc- 
tion of the flow (l6|. Inset: distorted zoom of a narrow strip 
containing the diagonal, magnified in the transverse direction 
to display the arc-shaped trajectory and steady state fluctu- 
ations. 

coupled by spin-conserving interdot tunneling with am- 
plitude t, exhibit an avoided level crossing as a function of 
detuning A. A uniform magnetic field B splits the (1, 1) 
triplet into the states | To ) , | T± ) with total z-projection 
of electron spin to = 0, ±1, respectively. 

The hypcrfine interaction gives rise to the Overhauser 
shift of the Zccman energy, which is different on the left 
and right dots: H^j^-^ = [g^isB + Asl{r))S^ . The 'po- 
larization gradient' As = sl—sr couples the states | Tq ) 
and I (1, 1)s ). Including this coupling, we find the energy 
levels £„ and eigenstates | n ) within the S^^^ = sub- 
space (spanned by the states | Tq ), | (1, 1)5 ), | (0, 2)5 )) 
by diagonalizing the Hamiltonian 

/ AAs/2 \ 

7^3X3 = ^As/2 t , (1) 

y t -~A~ihj/2 J 

where the energy of the unhybridizcd (0, 2) singlet state, 
£(o.2)s = ~A ~ *^7/2, includes an imaginary part that 



4 



accounts for its decay due to coupling to the continuum 
of states in the drain lead. Due to a large applied bias, 
we assume that the (1,1) state cannot decay directly to 
the continuum. After diagonalization, each of the states 
\n) obtains a nonzero | (0,2)5) component and decays 
with a rate 7„ = — 2Im[e„]/?i, see FiglSK. 

Current through the system results from electron 
transmission through any of the 5 states { | n ) , | T± ) } that 
may be populated when the second electron is injected 
from the lead. The total current 



(2) 



is determined by the inverse of the average of the lifetimes 
{r„.±} of these states, weighted by the probabilities Pn,± 
to load each of the states. 

Cotunneling or spin-exchange with the leads with rate 
Wcot adds an additional decay channel, leading to the in- 
verse lifctimes r-i = 7„ + W^cot, sec RefsU^lll-ii The 
inverse lifetimes of | T± ) are determined by the rates W± 
of resonant hyperfine flip-flop transitions to each of the 
states {\n)}, Eqs. ^ and (jH) below, and of cotunnel- 
ing: = W± + Wcot- We neglect spin-orbit coupling, 
which is weak in GaAs, but do not expect any qualitative 
changes if it is included. 

Figure [Sh shows current as a function of polarizations 
sl and sr on the two dots for a fixed set of external 
conditions. In the dark stripe along the main diagonal, 
Sl = sr, current is set by the cotunneling rate Wcot; 
away from this line the finite polarization gradient As 
mixes | Tq ) with the singlet states and gives rise to the 
large red plateau of enhanced current that spans most 
of the figure. The width of the dark stripe is set by 
the Overhauser field difference required to mix the | Tq ) 
state and the | (1, 1)5 )-like hybridized singlet state, and 
depends on tunnel coupling and detuning, as well as the 
cotunneling rate Wcot, which controls the saturation of 
current. Sharp yellow bands of further enhanced current 
appear where the Overhauser shift brings either | T+ ) or 
I T- ) into resonance with one of the states with S^^^ — 0. 
With the help of FigI3^, each line can be identified with 
a specific resonant transition. 

As seen in FiglSb, the {sl,sr) plane includes vast re- 
gions in which current is essentially independent of the 
value of polarization. In these regions, the mean value 
of polarization cannot be inferred and fluctuations about 
the mean do not induce fluctuations of current. In other 
regions, in particular near the line sl = sr, the deriva- 
tives dl/dsL^R are large; here current is very sensitive 
to small fluctuations of polarization which can result in 
large- amplitude intermittent fluctuations of current, as 
discussed above. 

As discussed in Reffioj. DNP arises during spin- 
blockaded transport as the result of competition between 
hyperfine decay of | T± ) , which changes the z-component 
of nuclear spin, and non-spin flip decay channels. Each 
time an electron in the state | !+(_) ) decays by hypcrflne 



spin flip, a nuclear spin is flipped from down to up (up 
to down). By adding the transition rates from \T±) to 
all three final states in the decaying S^^^ = subspace, 
obtained using Fermi's Golden Rule, we find the "bare" 
transition rates W±_i(jj) for hyperfine decay of | T± ) as- 
sisted by a nuclear spin fiip in the left (right) dot: 

_y^ |Qn±6„P A2(lTSfl)7« 

^'^ ^ 16Nr (e± - Rc£„)2 + (fi7n/2)2 ) 

with the hybridized states \n) = a„ | Tq ) -|-6„ | (1, 1)5 ) -f 
c„ I (0, 2)5 ) obtained from ([T]), and e± = g^isB ±A{sl + 
sr)/2. The states | n ) are characterized by a nonuniform 
electron spin density on the two dots, which introduces 
an asymmetry in the rates W-j- ^(/j) to flip nuclear spins 
on the left and right dots via the dependence on an, bn- 
The net spin flip rate in the left (right) dot is propor- 
tional to the total current /, Eq.Q, to the probability p± 
of loading \T±), and to the probability W±^l(^r)/{W±^l + 
W±^R + Wcot) that this state, when loaded, decays by hy- 
perflne spin flip in the left (right) dot, giving 



±,L{R) 



P±il/e) 



W±,L + W±.R + Wc. 



(5) 



Our separate treatment of | T± ) and | To ) is valid in 
nonzero field where the degeneracy of these states is 
lifted. Near zero field one should employ a 5 x 5 gen- 
eralization of Eq. ([1]) as in RefllTl. 

Finally, using Eq.([5]), and including relaxation with 
rate Froi, we arrive at the equations of motion for sl 
and sr: 



Sl 
Sr 



2(r+,L - 
2(r+,H 



T-.l)/Nl 

- T^^r)INr 



Frcl Sl 
Frcl Sr, 



(6) 



where sl{r) are defined as (A^+ — N^)/N in each dot. 



III. SIMULATIONS AND RESULTS 

Equation ^ defines a flow, illustrated by the arrows 
in FiglSt, that describes the smooth dynamics of mean 
polarization. However, polarization is actually stochasti- 
cally driven by a train of electrons passing through the 
system, and executes a directed random walk around the 
flow (|6]). Fluctuations arise from shot noise in the number 
of entering up and down spins, from the random compe- 
tition between hyperfine and cotunelling decay channels, 
and from nuclear spin diffusion/relaxation. We simu- 
late this random walk by stochastically loading electrons 
into each of the 5 transport channels and then generating 
a corresponding sequence of randomly distributed decay 
times and numbers of nuclear spin fiips, with mean values 
given by Eqs. (HI), dS]), and (g]). 

From the simulation we obtain current and polariza- 
tion trajectories as shown in FigslT}; andl^b, see appendix 
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ing from mixing of the triplet levels with the (1,1) singlet 
by the random hyperfine fieldii. However, when detun- 
ing is comparable to the tunnel coupling (dotted box), 
the traces show diminished zero field peaks flanked by 
noisy regions exhibiting large fiuctuations and stable re- 
gions of high current at higher fields. The boundary be- 
tween noisy and stable high current systematically moves 
to higher field as detuning is increased, and is hystcrctic 
in the sweep direction. 

By including a time-dependent external field in the 
simulation, we produced the field sweep traces shown in 
FigHb. The low-field boundary between noisy and quiet 
regions depends on detuning in a similar manner to that 
observed in the experiment, while on the high field side 
we find an additional current-step not observed within 
the range of available experimental data. Based on the 
corresponding behavior of DNP in the simulations, we 
thus interpret the transition from noisy to stable current 
in the experiment as an indication that the polarization 
quasi-fixed point, s^j^ = in Eq.(|ni), moves away from 
the sensitive region sl ~ sr- 



FIG. 4: Current during magnetic field sweep dB/dt < 0. a) 
Experimental field sweeps for several values of A, see Eq.([T]), 
in the regime described in Ref0, Fig. 3d (A = could only be 
approximately located experimentally). Curves are offset by 
0.1 pA for clarity, b) Simulated field sweeps in the field regime 
indicated by the dotted box in panel a), showing bistable 
current similar to experiment. In both cases, the threshold 
for noisy to quiet current moves to higher field with increasing 
detuning. To treat the case of B ~ 0, a 5 x 5 formulation akin 
to that of RefllTl must be used. 



for parameters. For dots of unequal sizes, 7^ Nj^, the 
fiow ^ is asymmetric with respect to and sr. In par- 
ticular, the rates ([3|) and (|4|) favor spin flips in the smaller 
of the two dots due to the increased hyperfine coupling 
per nuclear spin. As a result, the system can follow an 
arc-shaped trajectory like that shown in Figures [T]3 and 
ISt, in which polarization passes through the insensitive 
region sl ^ Sfi during its build up, eventually returning 
to a steady state sl ~ Sfj where polarization fluctuations 
result in large fluctuations of current. 

Slow fluctuations with a power spectrum close to 1//^, 
see Figl2]:, arc indicative of diffusion, which may be 
driven by nuclear dipole-dipolc interactions, as discussed 
in RcfH, or by current as described above. Unlike other 
sources of steady-state spin fluctuations, the shot-noise 
mechanism is intrinsic to spin-blockade and its inten- 
sity can be controlled by current. Alternative mech- 
anisms can thus be distinguished through the current- 
dependence of the underlying diffusion coefficient. 

Experimentally, current was also measured during slow 
sweeps of magnetic field (see Fig|4^). These data were 
obtained in the same regime as that of Fig. 3D in 
with a small change in tunnel coupling. At large A, cur- 
rent displays a simple peak at small magnetic fields aris- 



IV. CONCLUSIONS 

The mechanism described above, based on spin dy- 
namics driven by electron shot noise, provides a natural 
explanation for the systematic observation of regions of 
stable and strongly ffuctuating current. We propose that 
regions of high, stable current, see e.g. A = —lO^cV, 
B > 60 mT in FigHl indicate that the system tends to an 
asymmetric fixed point with a sizable difference between 
the hyperfine fields in the two dots. As recently demon- 
strated, such states can be used to perform controlled 
manipulations of the double dot electron spin states^,. 

We gratefully acknowledge financial support from 
FOM, NWO and the ERC (L.V.), NSF grants DMR- 
090647 and PHY-0646094 (M.R.), and the Intehigence 
Advanced Research Projects Activity (lARPA), through 
the Army Research Office (L.V. and M.R.). 



Appendix A: Simulation Parameters 

The main text describes the microscopic model used to 
generate the simulated current and polarization traces in 
Figures 1-4. The behavior exhibited by the model is sen- 
sitive to a number of parameters, many of which are not 
well characterized for the experimental system. While on 
the one hand the existence of many parameters makes di- 
rect comparison to experiment more difficult, it can also 
be seen as a necessary consequence of the fact that such a 
wide variety of complex phenomena have been observed 
in this system. To this end, we have attempted to in- 
clude the minimum number of ingredients necessary to 
produce the phenomena of slow quiet transients followed 
by steady state fiuctuations. We chose parameters with 
plausible values for realistic systems, see table HI which 
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led to simulated traces clearly demonstrating the phe- 
nomena of interest on field and time scales approximately 
comparable to those observed in experiments. 

In the table below, and Nb, are the numbers of 
nuclear spins in the left and right dots, Wcot is the co- 



tunneling rate, t is the tunnel coupling, A is the detuning, 
7 is the decay rate of (0, 2)s [see Eq.(l), main text], Fioi 
is the phenomenological relaxation rate for nuclear po- 
larization within the double dot, Eq.(6) of the main text, 
and B is the magnetic field strength. 
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TABLE I: Simulation Parameters 



Figure 


Nl 


Nr 


Wcot 


t 


A 


7 


Trcl 


B 


Ic 


4.8 X 10^ 


5 X 10^ 


2 X 10^ 




1 


6 ^eV 


3 fieV 


0.3 ^eV 


3.6 X 10"^ s"^ 


80 mT 


2b,c 


0.75 X 10'' 


1.25 X lO'' 


1 X 10* 




1 


2 ^eV 


4 /leV 


0.5 ^eV 


0.6 X 10"^ s-^ 


300 mT 


3c 


0.75 X 10^ 


1.25 X 10** 


1 X 10* 




1 


2 fieV 


4 fj.eV 


0.5 fieY 


0.6 X 10"^ s"^ 


400 mT 


3b" 


1 X 10^ 


2 X 10** 


1 X lO'' 




1 


15 fieY 


10 /xeV 


0.1 fieY 


N/A 


650 mT 


4 


1.6 X 10'^ 


2 X 10'' 


1 X 10^ 




1 


1.5 ^eV 


3-8 /xeV 


0.1 ^eV 


1.2 s"^ 


50 mT-350 mT 



"Parameters for Fig. 3b are chosen to most clearly display the 
sharp resonance features due to level crossings. 



